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ABSTRACT 

Within the internal shock scenario we consider different mechanisms of high energy 
(> 1 MeV) photon production inside a Gamma Ray Burst (GRB) fireball and derive 
the expected high energy photon spectra from individual GRBs during the prompt 
phase. The photon spectra of leptonic and hadronic origins are compared within dif- 
ferent sets of parameter regimes. Our results suggest that the high energy emission 
is dominated by the leptonic component if fraction of shock energy carried by elec- 
trons is not very small (e.g. e e > 10~ 3 ). For very small values of e e the hadronic 
emission component could be comparable to or even exceed the leptonic component 
in the GeV-TeV regime. However, in this case a much larger energy budget of the 
fireball is required to account for the same level of the observed sub-MeV spectrum. 
The fireballs are therefore extremely inefficient in radiation. For a canonical fireball 
bulk Lorentz factor (e.g. T — 400), emissions above ~ 10 GeV are attenuated by 
two-photon pair production processes. For a fireball with an even higher Lorentz fac- 
tor, the cutoff energy is higher, and emissions of 10 TeV - PeV due to 7r°-decay can 
also escape from the internal shocks. The flux level is however too low to be detected 
by current TeV detectors, and these photons also suffer attenuation by external soft 
photons. GLAST LAT can detect prompt emission of bright long GRBs above 100 
MeV. For short GRBs, the prompt emission can be only barely detected for nearby 
bright ones with relatively "long" durations (e.g. ~ 1 s). With the observed high 
energy spectrum alone, it appears that there is no clean picture to test the leptonic 
vs. hadronic origin of the gamma-rays. Such an issue may be however addressed by 
collecting both prompt and afterglow data. A moderate-to-high radiative efficiency 
would suggest a leptonic origin of high energy photons, while a GRB with an ex- 
tremely low radiative efficiency but an extended high energy emission component 
would be consistent with (but not a proof for) the hadronic origin. 
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1 INTRODUCTION 



The study of Gamma Ray Bursts (GRBs) has been one of the most interesting areas in astrophysics in the past few years. Ongoing 
observational and theoretical investigations are disclosing the physical origin, characteristics of these objects as well as bringing 
new puzzles to u s. EGRET detected high energy photons from five GRBs coincident with triggers from the BATSE instrument 
Jones et all 19961) . GRB 940217 wa s detected by EGRET independent of BAT SE trigger, which has extended emission and with 
the highest energy photon of 18GeV jHurlev et al J 1994n . I Gonzalez et al] Jiooj) discovered a distinct high energy component up to 
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200 MeV in GRB 941017 that has a different temporal evolution with respec t to the low e nergy component. Although even higher 
energy gamma rays/neutrinos have not been firmly detected from GRBs yet, Atkins et al. i 2000l) have provided tentative evidence 
of Te V emission from GRB 970417A. For a long time, GRBs have been identified as potential sources of ultrahigh energy cosmic 
rays i Waxmar] 1995 : Vietri|[l997 ). Within the standard fireball picture (e.g. Meszaros 2006), there are about a dozen mechanisms 
that can produce GeV-TeV gamma-rays from GRBs (e.g. Zhang 2007). More theoretical and observational efforts are needed to 
fully understand high energy emission from GRBs. From the theoretical aspect, it is essential to investigate the relative importance 
of various emission components to identify the dominant mechanisms under certain conditions. 

The high energy photon spectra expected from GRBs during the prompt and the afterglow phases have been derived 
by various groups. In the scenario of external shock model the high energy photon spectra during the early afterglow phase 
due to synchrotron and synchrotron se l f Compton (SSC) emis s ion by shock accelerated relativistic electron s and protons 
have been studie d (Mesza ros et al. 1994; Meszaros & Rees 1994 ; Panaitescu & Meszaros 1998 : Wei & LuJ 1998 : Totanil 1998 
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2007; Gou & Meszaros 2007). In the case of a strong reverse shock emission component, the SSC emission in the reverse shock 



region or the crossing inverse Co mpton processes between the forward and reverse shock regions are a lso important (Wan g et al 



Jw 

2001ijlbl: uVer & Waxmarj|2005h . The discovery of X-ray flares in early afterglows in the Swift era l IBurrows et al . 2005) also 
opens the possibility tha t scattering of the flaring photons from the external shocks can give strong GeV emission l Wang et"al] 
20061 : lFan&PirarJl2006h . The eff ect of cosmic in f rared background on high energy delayed 7-rays from GRBs has been also 
widely discussed i n the literature ( Dai & Lv] 2002 ; Steckei 2003 ; Wang et al. 2004 ; Razzaque et al. 2004 ; Casanova et al 2007 ; 
Muras e et al ] |2007h . The most important high energy emission component is believed to be emitted from the prompt phase. Swift 
early X-ray afterglow data suggest that the GRB prompt emission is of "internal" origin, unlike the external-origin afterglow emis- 
sion (Zhang et al. 2006, cf. Dermer 2007). The most widely discussed internal model of prompt emission is the internal shock model 
(Rees & M eszaros! 1994). Within the internal shock model the spectrum of high energy photons expected during the prompt phase 



has been studied i Pilla & Loeblll998l ; liragile et alj [20041 : iBhattachariee & Gupt<j|2003l ; iRazzaque et alj|2004l : IPe'er & WaxmarJ 
2004; Pe'er et alj|2006h . The various processes of high energy photon production in the internal shocks are electron synchrotron 
emission, SSC of electrons, synchrotron emission of protons, photon production through tv° decay produced in proton photon 
(P7) interactions and radiations by secondary positrons produced from n + decays. In this paper we consider all these processes 
self-consistently with a semi-analytical approach and study the relative importance of each component within the internal shock 
scenario. The derived photon spectra are corrected for internal optical depth for pair production, which is energy-dependent and also 
depends on various other parameters of GRBs e.g. their variability times, luminosities, the low energy photon spectra inside GRBs, 
and photon spectral break energies. If the electrons cool down by synchrotron and SSC emission to trans-relativistic energies, then 
they accumulate near a value of Lorentz factor of around unity. The accumulated electrons affect the high energy photon spectrum 
by direct-Compton scatterin g and other processes, w hich make the spectrum significantly different from the broken power laws 
considered in this work, see JPe'er et al.B2005ll2006l) for detailed discussions. In any case, for the values of parameters considered 



in the present p aper this effect is not s ignificant. 

GLAST's jGehrels & Michelsonll 1999b burst monitor (GBM) will detect photons in the energy range of lOkeV to 25MeV and 
large area telescope (LAT) will detect photons in the energy range of 20MeV and lOOOGeV. With a large field of view (> 2 sr 
for LAT), GLAST will detect high energy phot ons from m any GRBs and open a new era of studying GRBs in the high energy 
regime. This is supplemented by AGILE dLongo et al.ll2002h . which is designed to observe photons in the energy range of 10-40 
keV and 30MeV- 50GeV and also has a large fiel d of view. There are several other ground based detectors e.g. Whipple/VERITAS 
1 Horan et al.ll2007l) . Milagro dAfkins et al.l l2004). which have been searching or will search for ~ TeV photons from GRBs. De- 
tections or non-detections of high energy gamma rays from GRBs with space-based and ground-based detectors in the near future 
would make major steps in revealing the physical environment, bulk motion, mechanisms of particle acceleration and high energy 
photon production, photon densities, etc., of GRBs. 



2 ELECTRON SYNCHROTRON RADIATION 

We define three reference frames: (i) the comoving frame or the wind rest frame is the rest frame of the outflowing ejecta expanding 
with a Lorentz factor V with respect to the observer and the central engine; (ii) the source rest frame is attached to the GRB 
central engine at a redshift z; and (iii) the observer's frame is the reference frame of the observer on earth, which is related to the 
source rest frame by the redshift correction factor. We denote the quantities measured in the comoving frame with primes. The 
shock accelerated relativistic electrons lose energy by synchrotron radiation and SSC in the shock region. Assuming a power law 
distribution of fresh electrons accelerated from the internal shocks and considering a continuous injection of electrons during the 
propagation of the sh ocks, the relativi stic primary electron number distribution in the comoving frame can be expressed as a broken 
power law in energy 1 Sari et all 19981) 
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dN e (E' e ) J E' e v E' e m < E' e < E' e>c 



dE' e 1 E' e p E' e c < E' e 

in the case of slow cooling, where E' e m is the minimum injection energy of electrons and E' e c is the energy of an electron that loses 
its energy significantly during the dynamic time scale, known as the cooling energy of the electrons. If the electrons are cooling 
fast so that even the electrons with the minimum injection energy have cooled during the dynamical time scale, by considering 
continuous injection of electrons from the shock the comoving electron number distribution can be expressed as 

dN e (E' e ) f E'~ 2 E' e . c < E' e < K,m n , 

dE> K \ e'-*- 1 K, m < K ( ) 

If the electrons cool down to sub-relativistic energies then they accumul ate near electron Loren tz factor j' e ~ 1. This effect may 
distort the high energy photon spectrum by direct-Compton scattering dPe' er et al.ll2005i kOQfj) . and we focus on the parameter 



regime where this effect is not significant. The energies in the source rest frame and the comoving frame are related as E e ~ TE' e , 
where T is the average bulk Lorentz factor of the GRB fireball in the prompt phase. The expression for the minimum injection 
energy of electrons in the comoving frame is E' e m — m e c 2 -fpg(j>)^ s -^, where g(p) = for p >> 2 and g(p) ~ 1/6 



cj — — ~ o ~ —K,rrt iccvr/ me f£p > cr \r- / p _i r- - - - J \r / ~i~ 

for p — 2 i Razzaque & Zhang 2007), m p , m e are the masses of proton and electron, respectively, and 7 P m p c 2 is the average 



internal energy of protons in the comoving frame. We have assumed 7 P to be of the order of unity (in principle 7 P could be smaller 
than unity). The total internal energy is distributed among electrons, protons and the internal magnetic fields within the internal 
shocks. The fractions of the total energy carried by electrons, protons and internal magnetic fields are represented by e e , e p and es, 
respectively, where e e + e p + es = 1. We have assumed that all the electrons and protons are accelerated in internal shocks. In 
reality, the shock accelerated particles may be only a fraction of the total population and additional fractional parameters (£ e , £ p ) 
may be introduced (e.g. Bykov & Meszaros 1996). In such a case, the following treatments are still generally valid by re-defining 

e' e — e e /£ e and e' p = e p /£ p , while the relation e e + e p + es = 1 still holds. 

The relativistic electrons lose their energy by synchrotron radiation and inverse Compton scattering l Panaitescu & Meszaros! 



1998; Sari & Esin 2001; Zh ang & Mes zaros 2001). The comoving cooling break energy in the relativistic electron spectrum can be 



derived by comparing the cooling and the dynamical time scales. The comoving cooling time scale t' cool of electrons is a convolution 
of the cooling time scales for synchrotron radiation t' syn and for inverse Compton (IC) scattering t' IC 

111 

''cool 1 SS« L IC 

We denote U as the internal energy density of the internal shock, and U e , Ub as the energy densi ties of electrons an d magnetic 
fields, respectively. The energy density of the synchrotron radiation is U e ,syn = ?ly e = ^y 17 I Sari & Esinl 2001 ). where the 



radiation efficiency of electrons is rjc — [(-Ee.c/^e.m) P i 1] f° r sl° w an d fast cooling, respectively, and 



i+Y e i+y e 

2- 



Y _ Lc.IC _ Ue.syn _ -1 + yj 1 + 4?? e £ c JtB 
Le ;S yn Ub 2 

denotes the relative importance between the IC and the synchrotron emission component^ L e jc and L etSyn are the luminosities 
of radiations emitted in SSC and synchrotron emission of relativistic electrons respectively. The inverse of the cooling time scale 
of electrons can be expressed by the power divided by energy (E' e = m e 7eC 2 ), 

-, = -Oe,T/3 e 7e o(UB + U e ,sy n ) = -<T e ,T/3 e 7e J (1 + Y e ) , (5) 

t' cool 3 m e c^ 3 m e c 2 

where g^t is Thomson cross-section of electrons, /3 e ' ~ 1 is the dimensionless speed of the relativistic electrons. The comoving 
dynamical time scale is t' dyn ~ Tt v , where T is the average Lorentz factor of the GRB, and t v is the variability time in the source 
rest frame of the GRB, which denotes the variability time scale of the central engine. Throughout the paper, we assume that electron 
synchrotron radiation from the internal shocks is the mechanism that power the prompt gamma-ray emission in the sub-MeV band. 
However, for standard parameters within this scenario the cooling time scale of electrons is much shorter than the dynamical time 
scale of GRBs. As a result the flux density (^E JiS dN ~<^ i E ~< * s 3 . J below the coo ling break energy is proportional to E^J 2 and cannot 
explain the harder spectral indices observed in many GRBs fchisellini etallbood) . If the magnetic field created by internal shocks 



decays on a length scale much shorter then the comoving width of the plasma, then the resulting synchrotron radiation can explain 



1 Strictly speaking, such a treatment is valid for the IC process in the Thomson regime. However, this is also a reasonable approximation if the 
peak of the spectral energy distribution of the IC component is in the Thomson regime, which is generally the case for the calculations performed 
in this paper. 
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E ssa * 0.24 ^L^TrrT^B'^ = 0.69 ^L^-^T^ ^e B ,^ ( ^ ) ' (9) 



some of the broadband GRB spectra observed by Swift dPe'er & Zhang| |2006). In this case the effective dynamical time scale is 
shorter by a factor of f c than its actual value. Hence, the ratio of the cooling and the dynamical time scale can be expressed as 

— -fc (6) 

cool 

at the cooling energy E' e = E' e ^ c . The expression of the electron cooling energy in the comoving frame can be written as 

in' ' 2 2 3m e c 2 f c t v 2r!/c,2 

#e,c = 7e, c m- e c = m e c — — = 530keV- . (7) 

irt v a e ,TcUeB(l + Ye) Liso,5ies,-i(l + ie) 

Here and throughout the text the convention Qa; = Q/IO 11 is adopted in cgs units. In the above expression Li SO is the luminosity 
corresponding to the energy Ei ao carried by all particles and the magnetic fields in the shocks. It is a fraction of the wind (outflow) 
luminosity Li SO ~ r]L w , where r\ is the efficiency of converting the kinetic energy of the wind to the shock internal energy. 
The luminosity Lt so and internal energy U are related as U = Lj SO / ( A-KY 2 rj a 2 c), wh ere r; s = V 2 ct v is the internal shock 
radius. The synchrotron spectrum is a multi-segment broken power law l Sari et al. 1998b separated by several breaks, including 
the emission frequenc y from electrons with the m inimum injection energy, the cooling break frequency, and the synchrotron self- 
ab sorption frequency I Rvbicki & Lightmanll979h . In the internal shocks, the magnetic field in the comoving frame can be expressed 
as dzhang & Meszaroslbooj) 

B' * 4.4 x 10 B G« 1 6B,-i) 1/a L^, 81 rif 18 r a - 1 = 1.5 x 10 6 G ^E^if^l (8) 

1 2 tv,— 2 

where £ is the compression ratio, which is about 7 for strong shocks. The synchrotron self absorption energy (E ssa ) in internal 
shocks can be expressed as (Li & Song 2004; Fan et al. 2005; cf. Pe'er & Waxman 2004) 

2/7 

,1 + Ye, 

where L 7jS = Li 30 e e rj e / (1 + Y e ) is the isotropic gamma-ray luminosity due to synchrotron radiation. The cooling break energy 
-Ee jC and the minimum injection energy E' e ^ m of the electrons define two break energies in the synchrotron photon spectrum. The 
cooling break energy in the photon spectrum in the source rest frame is 

- r ?(%) 2£ ^ * !■« >< I*" 3 eVrJ = 2.8eVt„,- a . f ^ ^ (tW^ 

47rVm e c 2 / m e c 2 \i/»ao,5xes(l + Y e )J (ii SO ,5ies,-i) 3/2 V 1 + Y e J 

Notice that _E 7 , C very sensitively depends on V and some other parameters so that it could become a large value when parameters 
change. For example, for B' = 10 4 G, V = 400, f c = 500, L iso = 10 51 erg s" 1 , t v = 0.01s and e B = 0.1 we get E 1>c ~ 1.9 
MeV. The break energy in the photon spectrum due to the minimum electron injection energy is 

/L ^' U ''' 0.58MeVr 2 (^) 2 B^ = 8.5Mev(^) 2 te es ,-iL, so ,5i) 1/2 (r^,_ 2 )- 1 (11) 

Assuming E asa < -E 7 , TO ,s < E ltC , a the photon energy spectrum from synchrotron radiation of slow-cooling relativistic electrons 
is as follows 

... . j-, v ( E^ jS Essa <C -£^7,s ^ £/ 7 ,m,s 

2 diV 7 , s (£/ 7 , a j I 4 /3+{p _ 3)/ 2 ( p _3)/2 p p < „ n9 . 

T> s j 171 I J - J 'y,rn,s - L/ 7,s J -'-y,m,s \ J-''y,s _ -*-'"y,c,s \ az -/ 

dii T' s p4/3+(p-3)/2pl/2 „-(p-2)/2 „ ' p 

V. *- J, y,m,s J - / 7 ,e,s- C/ 'y,3 £-"y,c,s _ ^ 7 ,s 

In the case of slow-cooling electrons for very small values of e e (e.g. ~ 10 -3 , which is relevant when the hadronic emission 
component becomes important), the break in the photon spectrum due to the minimum injection energy of electrons goes below 
the synchrotron self absorption energy. The order in the spectral break energies becomes i5 7 . m , s < E SS a < E^ tC)S , and the 
spectrum is also modifi ed. The spectral indi ces of the electron synchrotron spectrum for different ordering of the spectral break 
energies are derived bv lGranot & Saril J2002b . For £ 7 , m , s < E 7 , s < E ssa the spectral index of E^ " 7 ^ 7 -'' is 7/2, and for 
E lt3 < E ltm ,s the spectral index is 3. The indices of the spectrum between E saa , -E 7lCjS and above -E 7 , c ,s remain as — (p — 3)/2 
and — (p — 2)/2, respectively. When E ssa is greater than both i5 7 , m , s and E ltC , s their relative ordering becomes unimportant. In 
that case the spectral indices of E 2 ^ dN ^^- E ~ l - s ' are 7/2 between i5 7im , s and E ssa , and — (p — 2)/2 above E S sa- Below J3 7im , s 
the index is 3. 

For fast-cooling electrons the synchrotron photon energy spectrum for E SS a < E ltC]S < -E 7 , mjS is 

. -l"y,s Essa < E JiS < E 7>C|S 

2 ClM 1 ,s(ti-/,s) I 5/6 pl/2 p p p n ~s 

p5/6 p(p-l)/2 p-(p-2)/2 p 
J- / 'y,c,s-*- / 'y ,m,s J-J^^s S2j<y 
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When the ordering of break energies in the photon spectrum becomes -B 7 , c ,s < E S3a < i5 7 , m , s the photon energy spectrum is 

J]\7" { 7~i \ f E "f> s Essa *C E"/,s ^ E^f^cs 

2 dN^, s {K 1>s ) I 9/g i/2 

7? s j T-i ^ \ J - /, y,c,s ■*- J, y,s -C / '7,c,s J- Jr y,s _ ^7 ,m,s V 1 ^"/ 

«^7,s p 9/8 p(f-l)/!p-(p-2)/2 p . p 

The total energy emitted in synchrotron radiation by relativistic electrons is Ei ao rj e e e /(l + Y e ). The normalisation constant for the 
synchrotron photon energy spectrum can be calculated from 

" B "' mQX riiV 7 , s (£ 7 , s ) _ r? e e e 

imiB d£ T , s d ^ s - £ » o (iTiy (15) 

The maximum electron energy E e , max can be calculated by equating the acceleration time and the shorter of the dynamical 
and cooling time scales of the relativistic electrons. The expression of the acceleration time scale is t' acc = 2nC,rL{E' e ) / c = 
2nC,E' e / eB' c. Here rL(E' e ) is the Larmor radius of an electron of energy E' e in a magnetic field B' , £ can be expressed as 
C ~ fish~ 2 y, where /3 s f, is the velo city of the shock in the com oving frame of the unshocked medium and y is the ratio of diffusion 
coefficient to the Bohm coefficient Rachen & Meszarosl 1 1998 ). In ultra-relativistic shocks f3 s h ~ 1 and numerical simulations for 



both parallel and oblique shocks gives £ ~ 1. With 

^acc ^^ n [^cool:^dyn] J (16) 

one can derive the maximum comoving electron energy 

K, ma * =rntn \s.5( —^^f^—) 1 ^ , 14.3 x 10 7 r 2 ^,_ 2 B^l GeV (17) 
L \Li a0}5 ie B (l + Y e )J J 

For electrons, the cooling term (first term in the bracket) always defines the maximum electron energy. The maximum synchrotron 
photon energy in the source rest frame can be then derived as 

„ _ r ^h f E' emax \ 2 eB' c „ , r / ^lB' 5 2 tl_ 2 \ „ / V 



. / E^\ eBc = Q 4gGeV / i2^_2 \ = 1Q2GeV / r2 \ 
4ttV m e c 2 J m e c 2 \L iao , 51 e B ,-i{l + Y e ) J \l + Y e J 

This is used in eqn.(|15b to define the normalization of the spectrum. The result has a very steep dependence on T. We also notice that 
B' is not an independent parameter, but can be calculated from other parameters according to eqn.®. For example, for V — 400, 
Liso = 10 51 erg/s, t v = 0.01s and es, e e ~ 0.1, the magnetic field is of the order of 10 4 G and the maximum photon energy 
becomes a few hundred GeV. 



3 ELECTRON INVERSE COMPTON SCATTERING 

The relativistic electrons can be inverse Compton scattered by low energy synchrotron photons inside the GRB fireball and transfer 
their energy to high energy photons. Below, we derive the IC photon spectrum using the electron and synchrotron photon spectra. 

dN^E^) _ 1 f dN^EJ f dN^E^) 

-lET dEe x J dEl , s dE ^ s (19) 

The electron Lorentz factor (je), IC and synchrotron photon energies (-E 7 ,i, Ej :S ) are related as E lt i ~ r y' e 2 E JtS , this can be used 
to simplify the above equation. The final expression for the IC photon spectrum considering slow cooling of electrons is 

E^ j E ssa ,i <C E^^i ^ ^? 7 ,m,i 



2 dN^(E^) s 4/3 + (p-3)/2 ( P -3)/2 
12 L \ ijyJ / 7,771,1 7,1 

J ^ dE^ K \ Et / t + } p - 3)/2 E i JLEZ<; p - 2)/2 l. , l l h 



^~t,i T7? K ) „4/3+fp-3)/2 ,-,1/2 „-(p-2)/2 „ , „ „ „ 



y,m,i 7,c,i 7,K 7,3 7;" 7i l 

Here J5 ssa ,i = ~f' e , m 2 Essa, E 1>m ,i = 'y'e,m 2 E-,, m ,s, and £ 7lC ,i = -y' e>c 2 E-, tC , a , where 7^ m = E' e}m /m e c 2 = 
g(p)(m p /me)(e e /ep), 7e, c = E' e ^/m e c 2 are Lorentz factors corresponding to the minimum injection energy of electrons and 
the cooling break energy of electrons. In the case of fast cooling E ltm ,i > E ltC ,i and the IC photon spectrum has to be modified 
accordingly. 

{E~ i E S sa,i *C E^^i ^ E~y fC ,i 

p 5/6 p p-l)/2 F -(f-2)/2 p p , p l. Zi J 

r;5/6 E (p-l)/2 E (p-2)/2 E -(p-2) „ „ 
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In eqn.(|2T) the expressions for E ssa ,i, E 1<Cti and Z? 7 , m ,j are E ssa ,i = j' e}C 2 E 3sa , E JlCli = 7e, c 2 -E 7 , c , s 
7e,m 2 ^7.m,s- When E e E lt3 » T 2 m 2 c 4 the cross section for IC scattering decreases as the scattering enters the Klein Nishina 
(KN) regime. A break in the photon spectrum at U 7] j = E^,k appears when the Klein Nishina effect becomes important. We define 
a parameter n = Y 2 m^lX' ' wnere E 1 , pea k = max[E-,,c,s; I? 7 , m , s ]. The KN regime starts when k = 1 (e.g. Fragile et al. 2004), 
and 

£ 7 ,k = = 2.5 GeV— (22) 

In the KN regime the emissivity of electrons decreases by k 2 , and the photon energy spectral index simply follows the electron 
energy spectral index, i.e. — (p — 2). The IC photon spectrum in eqn.J20b can be normalised as 

7> l 1 771 tA-C/ 7 ,i — £jiso 1 . -tr ) K^ J ) 

where E~,, ma ,x,i = TEL mnT due to the KN effect. 



4 PROTON SYNCHROTRON RADIATION 

Relativistic protons lose energy by synchrotron radiation and photo-pion (ir°, 7r + ) production inside GRBs. They interact with the 
low energy photons in the GRB environment and pions are produced. There is a threshold energy for this interaction (pr/) to happen, 
EpEj > 0.3GeV 2 r 2 , where E p and _E 7 are proton, photon energy in the source rest frame respectively. The 7r°s decay to a pair of 
high energy photons, while the 7r + s decay to neutrinos and leptons. The threshold condition therefore suggests that the photon-pion 
related high energy spectrum is typically more energetic than the electron IC spectrum. We assume that the proton spectrum in the 
internal shocks can be expressed as a power law in proton energy. We consider a proton spectral index similar to electrons for our 
present discussion. Since protons are poor emitters, we only consider the scenario of slow-cooling in the comoving proton spectrum 



dN p (E' p ) J E' p p E' p>m <E' p <E' P 



dEp 1 Ep p ^'p,c ^ 



p,c 



(24) 



where E' p m is the minimum injection energy of the protons and E' p c is break energy in the spectrum due to proton cooling. The 
minimum injection energy E' p m = r y' p m p c 2 g(p), where g(p) = 2^2 forp S> 2 and g{p) ~ 1/6 forp = 2. The cooling break 
energy can be derived by comparing the comoving and the cooling time scales. The inverse of the cooling time scale t' coo[ of a 
proton is 

1 1 1 

— = — + T (25) 

The photo-pion cooling t i me scale t'^ has been derived earlier in the context of estimation of neutrino fluxes from GRBs 



1 Waxman & Bahcalllll997l ; iGupta & Z hang 2007J). If f w is the fraction of proton energy going to pion production in the A res 



onance of jry interactions one has 1/^4 ~ f-w/t'dyn where the comoving time scale i£j t' dyn = Tt v . The peak value of p-f 
interaction cross section at the A resonance is o pl = 5 x 10~ 28 cm 2 . This is much higher than the Thomson cross section for 

protons <j Pt T = fe,T, where ct s ,t = 6.625 x 10~ 25 cm 2 . We therefore neglect the IC process of protons. Substituting for 

t' ayn and <4 in eqn.l!25t, we get 

VZrZ P ' TPp m P c 2 m P c* + TT v (26) 



where, f3' p is dimensionless speed of relativistic protons. We use the general expression for f w from lGupta & Zhang (2007) 



( 1.34°2-1 / E v xq 2 -l rp , Tp 

f( W ^_f) «2+l \E pb ) U p <U pb 



where 

^ _ 0.9-tnso,51 1 Vete ^g) 



2 In this definition, on average protons loose ~ 20% energy in the time scale oft' n . Although it is not strictly the e-folding timescale usually used 
to define cooling, for order-of-magnitude estimates this is good enough. 
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In our present discussion a 2 = (p + 2)/2 and a\ — (p + l)/2. E JtPea k t M e v is the peak energy in the electron synchrotron 
photon spectrum expressed in MeV, and Li so ,Bi is the GRB luminosity in unit of 10 51 erg s _1 , which is the typical value for GRB 
luminosities. E P b = 0.3r 2 /_B 7i p ea fc i G e ^GeV is the threshold proton energy for interaction with photons of energy -B 7 , pe afc,GeV- 
For typically observed values of GRB parameters one has E p b ~ 1 PeV. The break energy in the proton spectrum due to proton 
cooling can be calculated by comparing the comoving and cooling time scales of protons as discussed in the case of electrons in §2. 
We assume /3 P ~ 1 then for E p < E p b the expression of cooling break energy in the comoving frame is 

F' Js.fi 4 2cUeB \ k ^ a2 ~ 1 Y 1 _ 10 8 GeV/ c , 2 / L ISO , 51 6 B jp LW'^ y 1 

p ' c Tt v \3 p ' TPp mlc± + E pb Tt v a 2 + l J T 2 U,-2 \ T 2 H% t _ 2 E pb (J>eV)T 2 t v ,- 2 a 2 + 1 J ^' 

t' d 

where f c = ., — . The synchrotron photon spectrum from relativistic protons is 

cool 

2 dNy :PS (E^ ;Ps ^) ( E^^ ps ^ £'7,m,ps ^ E*y } p 8 -E'^.c.ps /in\ 

7 ' pa dE K I F 1/2 p-(P- 2 )/ 2 p <- p (iU) 

The minimu m injection energy in the p hoton spectrum from proton synchrotron radiation is related to that from electron synchrotron 
radiation as dZhang & Mesz aros 2001) 

E I E' \ 2 / ro \ 3 

s - in f,m,ps / p,m \ / " L e \ (31) 

£^7,m,s V E e m ) V 7Tlp ) 

The cooling break energy in the photon spectrum from proton synchrotron radiation is the characteristic synchrotron photon energy 
for proton energy E' p c . To normalize the proton synchrotron spectrum, it is important to find out the relative importance between 
proton synchrotron radiation and p7 interactions. Similar to the treatment of electrons, one can define 

Y _ -^P.P7 _ °P7 Ue,syn _ °"p7 y ^2) 
Lp : syn a p,T Ub a P-T 

^PyP'y and L p s y n 

tively. Notice that protons interact with the synchrotron emission of the electrons, so that V e enters the problem. Eqn. J32b suggests 
that Y p is usually much greater than unity since a pl ^> <j v ,t- As a result, most of the proton energy is lost through p7 interaction 
rather than proton synchrotron radiation. 

The proton synchrotron photon spectrum can be normalised as 

„ dN~ llPS {E ltPS ) JJ7 , _ p tpflp 

Qi£-j'y > ps -L ~r J- P 



where rj p — \E' pc / E' p ^ m j . The maximum proton synchrotron photon energy is derived by E ltmax ,ps = 

^ E y Tl °T ^ ^ '-p , where E' p ^ max is again defined by comparing the comoving acceleration time with the shorter of the co- 
moving dynamical and cooling times scales 

D/p6j.2 . 1/2 

iiso,5ies(l + Y p ) , 



E 



50 (7 nl vJ ,1.4x 10 6 r 2 t„,- 2 B6 



TeV . (34) 



E 



ft \^ Vlt v ,. 2 \^ 208 4 1/2 



TeV . (35) 



5 7T DECAY 

The relativistic protons interact with the low energy photons and photo-pions (7r°,7r + ) are produced as a result. The probabilities 
of 7T° and 7r + production are 1/3 and 2/3, respectively. Pions subsequently decay, i.e. it — > 77 and 7r + — > — > u fl i>^u e e + . 
As the cross section for the 77 interactions is much higher than the peak value of p7 interaction cross section, above the threshold 
energy of pair production 77 interactions are expected to dominate over P7 interactions. If the photon energy is 2m e c 2 ~ 1 MeV in 
the comoving frame, then in the source rest frame it is of the order of a few hundred MeV as the Lorentz factors are typically of the 
order of few hundred for canonical GRBs. For example, for T = 400 the photons of energy 400 MeV can produce photo-pions by 
interaction with protons of minimum energy E p ~ 120 TeV. The tt typically carries 20% of the proton's energy and the photons 
produced in it decay share its energy equally. Hence, the minimum energy of the photons produced from tv° decay is expected to 
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be ~ 10%-Ep ~ 12 TeV. The photon spectrum produced from ir° decay has been derived below using the proton spectrum defined 
in eqn.d24t and assuming the fraction / w /3 of protons' energy goes to ty°s. 

where, E~ v o c = 0.1_B P , C . For the expression for /„■, see eqn,l !27b . which contains a break energy. The break energy in the photon 
spectrum contained within f n is £L n o i = 0.03T 2 /ebr,GeV GeV assuming 10% of the proton's energy goes to the photon produced 
via 7T° decay, ebr is the break energy in the low energy photon spectrum (in the scenario of slowly cooling electrons it is the cooling 
break energy in the photon spectrum and for fast cooling electrons it is the photon energy corresponding to the minimum injection 
energy of electrons). The photon flux can be normalised in the following way 



E 



dN lw o(E 1 . K o) _E iao t p rj p Y p 



dE~o dE ^° - — T+Yp (37) 

where E„ ^o >min = 30F GeV and E 1 >5r o >max = Q.lE Pim ax- Although high energy photons (~ TeV) are absorbed by lower 
energy photons and e + e~ pairs are produced, at extreme energies the pair production cross section decreases with increasing energy 
i Razzaque et al . 2004). ] lence, ultrahigh energy photons can escape from the internal shocks for suitable parameters depending on 
the values of their various parameters and the low energy photon spectra. 



6 SYNCHROTRON RADIATION OF POSITRONS PRODUCED IN vr+ DECAY 

The shock accelerated protons may interact with the low energy photons to produce 7r + s along with 7r°s as discussed in the 
previous section. The 7r + s subsequently decay to muons and neutrinos. The energetic muons decay to positrons and neutrinos 
(P7 — > 7r + — ► — * e + UfiV^Ve). The charged pions, muons and the positrons are expected to lose energy through synchrotron 
radiation and IC inside the shock region. As the Thomson cross section for positrons is much larger than pions or muons, they 
are expected to emit much more radiation compared to the heavier charged particles. On the other hand, since these positrons are 
very energetic, most IC processes happen in the Klein Nishina regime. We therefore neglect the contribution of the positron IC 
processes. The positron synchrotron spectrum produced in pj interactions can be derived in the following way. The fraction of the 
protons' energy tranferred to pions is denoted by f n (eqn |27| ). If we assume that the final state leptons share the pion's energy 
equally then one fourth of the pion's energy goes to the positron. The energy of the positron spectrum dN d g°^ - at the energy E e + 
can be expressed using the proton spectrum defined in eqn. J24t 

, 2 dN{E, A ) 2 /,(£.,) I E'-' E„ < E„ 



where, E e + c is the cooling break energy in the positron specrum and 

ai + 1 (Wj;> h <=+ > ^e + b 

where /o has been defined in eqn.{|28l>. E e + b = 0.05E p b, E p b = 0.3 GeVF 2 /e(, ri GeV, and ebr,GeV is the break energy in the 
photon spectrum as defined earlier. The positron spectrum in eqn.l|38t can be normalised using the total energy carried by the 
positrons, 

dN(E e+ )(E e+ ) ^ 1 e pVp Y p E lao 
dE e + 6 1 + Y p 

The maximum and minimum positron energies are E e + max = 0.05E P: max and E e + min = 15F GeV (which is ~ 6TeV for 
r = 400). The synchrotron photon spectrum from the positrons can be subsequently derived using the same treatment for primary 
electrons as discussed in §2. The IC emission is in the KN regime and therefore not important. Also, photons having energies 
above a few hundred GeV are annihilated by lower energy photons as discussed in the following section. The relativistic muons 
produced in tt + decay lose energy by synchrotron radiation. We compare the decay and synchrotron energy loss time scales of the 
high energy muons. The maximum energies of positrons can be calculated in this way. If the muons decay before losing energy 
significantly high energy positrons are produced carrying approximately 5% of the initial proton's energy. On the otherhand if 
the muons lose energy before they decay lower energy positrons are produced. These positrons radiate energy and produce lower 
energy photons. The muons initially carry approximately 10% of the relativistic protons' energy hence, we expect the low energy 
photon flux produced by cooling of positrons is lower than that produced by relativistic electrons if e e and e p are comparable. 
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7 INTERNAL PAIR-PRODUCTION OPTICAL DEPTHS OF HIGH ENERGY PHOTONS 

Inside GRBs high energy photons interact with low energy photons to produce electron-positron pairs (e.g. Baring & Harding 1997; 
Lithwick & Sari 2001). The optical depth d epends on the values of various parameters of the GRB fireball. We follow the approach 
discussed in Bhattachariee & Guptal 1 2003) to derive internal optical depths of GRBs in detail. For two photons (a high energy 



photon 7^ and a low energy photon 7;), the pair production cross section depends on the energies of the photons and the angle 
between their directions of propagation. The cross section is 1 Berestetskii et al ]fl982 ) 



a lh r ri {E-y h ' , 0) 



3 n 







(3-/3' ) In 



1 + /3' 
1-/3' 



2/3'(2-/3' 2 ) 



(41) 



where ot is the Thomson cross section, and /3' = [1 — (£?' ^h/E^ )] 1 ' 2 is the center of mass dimensionless speed of the pair 
produced. The threshold energy of pair production with a high energy photon of energy Ely is 



E' lu th - 



2(m e c 2 ) 2 
E!y h (l - CO, 



(42) 



For the photons with energy higher than the threshold energy, the pair production cross section decreases with increasing photon 
energy 1 Jauch & Rohrlichl 1955 ; Razzaque et aljEo04 ). In the present work we calculate internal optical depths in different energy 
regimes using the cross sections with different energy dependences. The mean free path for jh 7; interactions l lh ll can be calculated 
using the low energy photon spectrum. 



r 1 



and, 



dE' 



E' 

11 • 



dn ll (E'-y, 

dEL, 



0) 



(43) 



d(cos8)(l-cos9)r' e (E' Jh ,8) 



(44) 



where 



is the specific number density of low energy photons inside the GRB. The low energy photon spectrum is ob- 



servationally known, as revealed by gamma-ray detectors such as BATSE and Swift. Theoretically, it corresponds to the electron 
synchrotron component as discussed in §2, which is a broken power law spectrum separated by the synchrotron self absorption 
break, the minimum injection break and the cooling break. The low energy photon flux is related to the observed luminosity 
through 



E~, 



, dn ll (E' ll ) 
dEL, 



dE~y, — TJsy — 



A-Kcr is 2 T 2 



(45) 



where L lt i SO is the isotropic 7-ray luminosity. We have taken it to be equal to the luminosity of the synchrotron photons emitted by 
electrons: L-y,i SO = L e , S yn = 6e ^f y p ■ In eqn.d44b we have three variables: ang le and photon energies E'„ n , E' lh . To simplify 
the integration in eqn.J44t we transform the integral with a new variable following lGould & Schreden ( 119671) 



E ll 



= s e 



2(m e c 2 ) 2 

with s = f"" E 2& , 
of the new variable s. It is then possible to write eqn.J44l( as 



(46) 



and G = h(l — cosO). As /3' = (1 — 1/s) 1 ^ 2 , the pair production cross section can be expressed as a function 



( mjc 
V EL, 



E-, 



idny^E'^) 



dE' 



dE: 



[so(E' lt )] 



(47) 



where 



?[so(K)l = 



r(s)ds . 



(48) 



and a(s) 



■ "inn 



For moderate values of s we use a(s) ~ 1 and for s » 1 it can be approximated as er(s) ~ ln(s)/s. 



The expressions for Q[sq(EL )] are (s§ — l)/2 and so(lnso — 1), respectively, in the two cases. Substituting for Q[so(E! Jl )] in 
eqnJ47t we derive the final expression for l^ h ~ fl (EL ). The internal optical depth Ti n t{E' lh ) is the ratio of comoving time scale 
and the mean time between two pair production interactions. 
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nnt(E' Jh ) = r ^cl-^(E', h ) (49) 

The final photon energy spectrum to be observed on Earth from nearby GRBs (neglecting further attenuations with the infrared 
background and cosmic microwave background) can be obtained by correcting the original flux for the internal optical depth and 
the redshift z of the source 

dA ^ 7 , o6 (g 7 , o6 ) = 1 E 2 dNj(E j) 
dE lt0b ~A-k<P z (1 + z) 1 dEj 

where 



E -t'°* TS^, - i ^\ E i— TB exp(-T m t(S 7 )) , (50) 



dz = JL f d *' (51) 

is the comoving distance of the source, Ho = 71km s _1 Mpc -1 is the Hubble constant, and Qa = 0.73 and Q m — 0.27 are 
adopted in our calculations. 



8 PHOTON SPECTRUM FROM SECONDARY ELECTRONS AND POSITRONS 



The secondary pairs carry a significant fraction of energy in the primary spectrum, and this energy is re-radiated and converted 
to photons. A more realistic treatment sh ould consider a photon-pair cascade process, which requires numerical calculations 
i Pe'er & Waxmanlliool ; be'er et alJlioOol) . iere instead we estimate the emission from the secondary pairs. We first calculate 



the photon energy spectra generated by different physical processes as discussed earlier. The photon spectra are then corrected 
for internal optical depths and subsequently the total energies carried by these photons are calculated by integrating the corrected 
photon energy spectra over photon energies. If we subtract the total energies carried by these high energy photons from their intial 
energies before including the effects of internal optical depths, we get the energies of the secondary e~ and e + produced in 77 in- 
teractions. These pairs are expected to have spectral indices similar to the high energy photons. With the knowledge of their spectral 
indices and the total energies carried by them the synchrotron photon spectra radiated by these secondary leptons are calculated. 
For the parameters adopted in this paper, it turns out that the emission contribution from the secondaries is below the emission level 
of the primaries, and hence, does not significantly modify the observed the spectrum. We therefore do not include this component 
in Figs. 1-5 , but caution that such a feed back process could be potentially important for the parameter regimes with high opacity. 
We refer to Pe'er & Waxmar] 1 2004 ) and Pe'er et al. ( 2006) for more detailed treatments of such cases. 



9 SYNTHESIZED SPECTRA AND DETECTABILITY 

Using the procedure delineated above, we have calculated the broad-band emission spectrum from internal shocks for a wide range 
of parameter regimes. In particular we focus on the various high energy emission components discussed above and their relative 
significance. Our results are presented in Fig. 1-5. In each set of calculations we have presented the internal optical depth after the 
final photon energy spectru m. For particles accelerated by ultra-relativistic shocks the spectral index is expected to be about 2.26 
Lemoine & Pelle tier ( 2003). Afterglow modeling suggests a larger scatter of p values for relativistic shocks, but p = 2.3 is close the 
mean value of the data (Panaitescu & Kumar 2002). In all our calculations, the spectral indices of relativistic electrons and protons 
are both assumed as p — 2.3. 

Figures 1-4 are the calculations for a typical long GRB with duration T90 = 20 s at redshift z — l(10sin the source 
rest frame). Since we do not know the physical condition of the internal shocks from the first principles, we vary the parameter 
regime in a wide range. In each set of calculations, we design the parameters to make the electron synchrotron emission peaking 
at the sub-MeV range (~ 0.36 MeV, 0.13 MeV, 0.6 MeV, 0.25 MeV for Figs. 1-4, respectively), as suggested by the data. The 
global energetics of the GRB is also adjusted so that the gamma-ray luminosity in the sub-MeV range is about 10 51 ergs s _1 as 
suggested by the observations. The variability time scale for these calculations is taken as t v = 0.01 s. The bulk Lorentz factor 
is adopted as T = 400 in Fig. 1-3, as suggested by the recent early optical afterglow observations (Molinari et al. 2006). In order 
to check how T affects the spectra, we also calculate the case of T = 1000 for the parameter set of Fig. 1, which is presented in 
Fig.4. In all the figures, the different components of the photon energy spectrum from a GRB for both electrons (e) and protons 
(p) are displayed with different line styles/colours. The observed energy fluxes E'l ob dN 7M E 7 ,ob) in umt of ergs/cm 2 

sec are 

plotted against the observed photon energy E Jt0 t(eV). The green long dashed curves represent the synchrotron emission from the 
relativistic electrons. The short dashed curves (blue) represent the IC spectrum from energetic electrons; the dash-dotted curves 
(light blue) represents the synchrotron emission of the relativistic protons; the triple short dashed curves (orange) represent for the 
synchrotron emission of the relativistic positrons produced in ty + decays; the ultrahigh energy emission component from ty° decays 
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is shown by the double short dashed curves (black) in the extremely high energy regime. The thin black solid lines represent the 
synthesized spectra of various components without including the effect of pair production attenuation. Depending on parameters, 
the pair opacity becomes important in the GeV - TeV range. The thick black solid lines represent the final photon spectrum after 
including the internal optical depths. In order to check whether the predicted high energy components are detectable by GLAST, 
we also plot an indicative GLAST sensitivity threshold in the 100 MeV - 100 GeV energy range. The GLAST sensitivity estimate 
is based on the criterion of detecting at least a few photons in the band based on the average effective area and photon incoming 
zenith angle of LAT. Background is negligible for GRB detections. This gives a rough fluence threshold of ~ 2 x 10~ 7 erg cm~ 2 
(B. Dingus, 2007, personal communication). The flux thresholds adopted in all the figures are therefore derived from the observed 
durations. For Tgo = 20 s, this gives a flux threshold of ~ 10~ 8 erg cm -2 s _1 . The sensitivity of VERITAS to photon above 
energy 200GeV has also been shown in our figures with pink dotted line. It is 2 x 10 _8 erg cm -2 s -1 (D. Horan, 2007, personal 
communication). 

Figure 1 is a standard "slow-cooling" leptonic-dominant case. The shock equipartition parameters are e e = 0.4 and es = 0.2. 
The isotropic shock luminosity is Li SO — 10 52 erg s _1 . The slow cooling factor f c = 2500 is adopted, which suggests that the 
post-shock magnetic field decays on a length scale shorter than the comoving scale (Pe'er & Zhang 2006). The thick black line 
shown on the right side around 10 15 eV is the it component after including the effect of absorption due to pair production, indicating 
the reduction of pair opacity at high energies (Fig. lb, see also Razzaque et al. 2004). In this figure the break energies in the photon 
energy spectrum appear in the order of E 3sa < £ 7 , m < E^, c in the electron synchrotron and IC spectral components. The spectral 
index of the photon energy spectrum is 4/3 between E ssa and i5 7 , m , — (p — 3)/2 between E lt , n and U 7]C , and — (p — 2) /2 above 
E JiC . Since e e is large, the leptonic components are many orders of magnitude stronger than the hadronic components. The value 
of Y p is much larger than 1 , so that the proton synchrotron component is below the components due to tt° decay and positron 
synchrotron radiation. 

We vary the values of the equipartition parameters (e e , e_e, e p ) and study the variations in the photon energy fluxes generated 
by various processes. The emission level of the electron IC spectral component decreases with decreasing e e (fixing e_g) since F e 
is decreasing. Moreover, as we decrease e e the minimum injection energy of electrons -E 7 , m also decreases. In the slow cooling 
regimes, it is E ltC that defines the peak energy in the electron synchrotron spectrum, which could be adjusted to the sub-MeV range 
by adopting a suitable f c value. The change of E Jim therefore mainly affects the calculated internal optical depth. 

By lowering e e , we check the parameter regime where the hadronic component becomes comparable. Since eletrons are much 
more efficient emitters than protons, the parameter regime for the hadronic component to be comparable to the leptonic component 
in the high energy regime is e e /e p ~ m e /m p < 10~ 3 H A similar conclusion has been drawn for the external shocks (Zhang & 
Meszaros 2001). In Fig.2, with e e = 10~ 3 , e s = 0.05 and e p = 0.849. In order to adjust E 1 , c to the sub-MeV range, f c = 50000 
is needed. In order to match the observed MeV emission flux by electron synchrotron, a large energy budget i s needed due to a 
small e e : Ei SO = 10 56 ergs and Li SO = 10 55 erg s _1 . Such a large energy budget has been suggest ed before ^ Totani 1998 ), bu t 
afterglow observations and modeling in the pre-Swift era have generally disfavored such a possibility l Panaitescu & Kumarll2002 ). 
In the Swift era, however, a large afterglow kinetic energy for some GRBs is not ruled out. For example, the bright afterglow 
of GRB 061 007 de mands a huge kinetic energy if the afterglow is produced by isotropic external shocks JMundell et al.l 120071 : 
Schadvetal1l2007l) . Modeling some X-ray afterglows below the cooling frequency requires a low es and/or a large afterglow 
kinetic energy at least for some GRBs Jzhang et alj2007l) . We therefore still consider such a possibility. In Fig.2, the break energy 
in the photon energy spectrum due to the minimum injection energy of electrons is below the synchrotron self absorption energy. 
The break energies appear in the order of E Jim < E 3aa < -E 7 , c in the synchrotron and IC electron spectra. The spectral index of 
the photon energy flux is 7/2 between E^, m and E asa , — (p — 3)/2 between E S3a and E ltC , and — (p — 2) /2 above E ltC . We can 
see that in the TeV energy regime beyond the maximum electron synchrotron energy, the positron synchrotron emission from 7r + 
decay becomes dominant. Moreover, when e e is small, Y e is small, hence Y p becomes small. In this case the proton synchrotron 
component becomes comparable to the spectral components due to synchrotron radiation of the secondary positrons and ir° decays. 
The internal optical depth is plotted in Fig. 2b, which peaks at a higher energy than that in Fig. lb. 

If the post shock magnetic field does not decay within a short distance (f c = 1), internal shocks are in the standard fast-cooling 
regime. We calculate such a case in Fig.3. The shock parameters are e e = 0.6, es = 0.2, Li SO — 10 52 erg s _1 , Ei so — 10 53 erg. 
In this case the break energies appear as in the order of Ec < E ssa < E m . The photon energy spectral indices are 13/8, 1/2 and 
— (p — 2)/2, respectively, in the three energy regimes. 

The pair opacity depends on the bulk Lorentz factor. When T is large enough, the ultra-high energy photons would have lower 



3 Pro ton energy loss and their contribution to high energy photon emission in the early afterglow phase has been studied earlier bv lPe'er & Waxmar] 
1 2005). Our results for the prompt emission phase are generally consistent with them. In order for the proton synchrotron component to be significant, 
even smaller e e (than 10 — 3 ) is demanded. Considering that photon-pion emission is more efficient than proton synchrotron emission, the condition 
e e /e p ~ m e /m p < 10 can allow the hadronic components to be comparable to (but not dominant over) the leptonic components. 
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Figure 1. A leptonic-component-dominated slow cooling spectrum, (a) The different components of the photon energy spectrum from the internal 



shocks for the following parameters in the slow-cooling regime: Ei, 



10"erg, L, 



10 52 erg/s, tv = 0.01s and f c = 2500. The thick solid 



black curve represents the final spectrum after including the effect of internal optical depths. The thin solid black curve represents the synthesized 
spectrum before including the effect of internal opical depths. The long dashed (green) curve is the electron synchrotron component; the short dashed 
(blue) curve is the electron IC component; the double short dashed (black) curve on the right side is for 7r° decay component; the triple short dashed 
(orange) line represents the synchrotron radiation produced by positrons generated in 7r+ decays; the dash-dotted (light blue) line represents the 
proton synchrotron component. The tiny red horizontal line between 10 s and 10 11 eV represents GLAST's threshold. The pink dotted horizontal 
line above 2 X 10 11 eV represents the sensitivity of VERITAS experiment (b) Internal optical depths plotted against energy for the parameters 
adopted in (a). 



internal optical depth and may escape from the internal shocks dRazzaque et al . 2004). To test this, in Fig.4, we re-calculate with the 
parameter set for Fig. 1, but increase T to 1000. The slow-cooling parameter f c is adjusted to 50 to maintain the sub-MeV energy 
peak. The results indeed suggest that the attenuation of the high energy photons is weaker. 

The observational breakthough in 2005 suggests that at least some short GRBs are low-fluence, nearby events that have a 
distinct progenitor than long GRBs (Gehrels et al. 2005; Bloom et al. 2006; Fox et al. 2005; Villasenor et al. 2005; Barthelmy et al. 
2005; Berger et al. 2005). To check the prospect of detecting short GRB prompt emission with high energy detectors such as GLAST, 
we perform a calculation for the parameters of a short GRB in Fig.5. Due to their short durations, short GRB detections are favorable 
for high luminosity and relatively "long durations". We therefore take an optimistic set of parameters with Li SO — 10 51 erg s - , 
Tgo = 1 s, and z = 0.1. Other parameters include: T — 800, t v = 1 ms, e e = 0.4, e_g = 0.2, e p — 0.4, f c — 50. The photon 
flux from synchrotron radiation of electrons peaks at O.lMeV. Fig. 5a suggests that the high energy component of such a burst is 
barely detectable by GLAST. The internal optical depth of this set of parameters does not grow to very large values (maximum 10), 
so that the attenuation signature is not significant in Fig. 5a. The dip around several 10 13 eV corresponds to the optical depth peak, 
above which the attenuated flux starts to rise. The abrupt drop at several 10 14 eV corresponds to the disappearance of the electron 
IC component at high energies. 
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10 CONCLUSIONS AND DISCUSSION 



We have calculated the broad-band spectrum of GRBs from internal shocks for a wide range of parameter regimes. We did not 
take into account the external attenuation of TeV photons by the infrared radiation background and that of the PeV photons by the 
cosmic microwave background. These external processe s would further attenuate our ca l culated spectrum in high energy regimes , 



and reprocess the energy to del ayed diffuse emission lDai&Lull2002l ; Isteckejliool ; IWang et al. 2004; Ra zzaque et al 



2004 



Casanova et al 2007; Murase et al. 1 20071) . Such processes are not relevant for most of the calculations presented, however, since the 
internal attenuation already cuts the observed spectrum below TeV. They are however important for high Lorentz factor cases in 
which more high energy photons are leaked out of the internal shock region. The external attenuation is also prominant for high 
energy emission from the external reverse/forward shocks and the external IC processes related to X-ray flares. These processes 
have been extensively discussed in other papers (referenced in Introduction) and they are not discussed in this paper. For nearby 
GRBs (e.g. z < 0.3), TeV emission is transparent. It is possible that ground-based Cherenkov detectors such as VERITAS, Milagro 

would detect TeV gamma-rays from nearby energetic GRBs. 

In previous treatments of hadronic components from internal shocks l Fragile et al. 2004 ; Bhattachariee & G upta 200^), the 
shock accelerated protons are assumed to carry m p /m e times more energy than electrons. This effectively fixed e c ~ m e /m p , 
which is not justified from the first principle. In this paper we have taken all the equipartition parameters e e , e p and es as free 
parameters, and explore the relative importance of various components in different parameter regimes. The dominant hadronic 
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Figure 3. A leptonic-component-dominated fast-cooling spectrum, (a) The spectra of various components. Parameters: e e = 0.6, eg = 0.2, 
e p = 0.2, tv = 0.01s, f c = 1, E iBO = 10 53 erg and L iso = 10 52 erg/s. Same line styles have been used as in Fig. 1 . (b) The corresponding 
internal optical depths. 



component emission becomes interesting only when e e is extremely small. Given the same observed level of sub-MeV spectrum, 
the total energy budget of the GRB needs to be very large. Inspecting the calculated spectra for different parameter sets (Figs. 1-4), 
one finds that there is no clean picture to test the leptonic vs. hadronic origin of the gamma-rays. Such an issue may be however 
addressed by collecting both prompt and afterglow data. A moderate-to-high radiative efficiency would suggest a leptonic origin 
of high energy photons, while a GRB with an extremely low radiative efficiency but an extended high energy emission component 
would be consistent with (but not a proof for) the hadronic origin. 

T he prompt emissio n produced by leptons including the effect of pair production has been discussed by Pe'er & Wnxmanl 

( 2004); Pe 'eretalJ ^2006). They calculated the emergent photon spectra for GRBs located at z = 1. The lower cut-off energy in 
the photon flux produced by leptons is determined by the synchrotron self absorption energy, the minimum injection energy or the 
cooling energy depending on the values of the various GRB parameters. Our leptonic-component-dominated cases are consistent 
with their results, although we do not explore cases with very high compactness. If the electrons cool down to trans-realtivistic 
energies then their high energy spectrum significantly deviat es from broken power law dPe'er et alj |2006, 2005). For our choice 
of values of the GRB parameters this effect is not important. iRazzaque et alj ( 120041) estimated the internal optical depth for pair 
production and showed that at PeV energies the optical depth decreases with increa sing photon energies. We have rederived the 
optical depths for values of GRB parameters. The results are generally consistent with dRazzaque et al .l2004h except that the growth 
of optical depth with increasing energy is more gradual before the optical depth peak. This is a result of including the whole low 
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Figure 4. The case of a higher Lorentz factor, (a) The spectra of various components. Parameters: T = 1000 and f c 
are the same as in Fig. 1. (b) The corresponding internal optical depths. 



50. All the other parameters 



energy photon spectrum (rather than the threshold energy photons) for calculating the pair production optical depth. The optical 
depths depend on the cross section of 77 interactions, the low energy photon spectra, the various break energies in those spectra, 
luminosities, variability times and the GRB Lorentz factors. A change in values of any of these parameters may affect the values 
of the optical depths at various energies. For high bulk Lorentz factors, the 7r° component may appear in the final spectra due to 
the reduced optical depths around PeV energ y. However, th ese ultra-high energy photons will be immediately absorbed in the GRB 
neighborhood by cosmic microwave photons (Stecker 2003). The reradiated energy by the e + e~ pairs would nonetheless contribute 
to the diffuse high energy 7-ray background (Cas anova et a]|2007l) . 

Upcoming 7-ray detectors have a good chance of detecting prompt emission from GRBs and reveal their physical nature 
during the prompt phase. Detection of the hadronic components is difficult but it would be possible to infer the dominance of these 
components by a coordinated broadband observational campaign if they are indeed important. More generally, detection or non 
detection of high energy photons in the prompt phase would constrain the values of various GRB parameters. In particular, the pair 
attenuation feature would help to constrain the bulk Lorentz factor of the fireball. Compared with EGRET, GLAST has a 10 times 
larger collecting area and a larger field of view. It is expected that GLAST LAT would detect high energy emission from a large 
number of bursts (mostly long GRBs and some bright, relatively "long" short GRBs), which will open a new era of studying GRBs 
in the GeV-TeV regime. On the other hand, it is difficult for VERITAS to detect prompt high energy gamma-rays even under the 
most optimistic conditions. High energy emissions from the external shock at the early afterglow phase for nearby GRBs may be 



© 2007 RAS, MNRAS OOO.Q]-?? 



16 Nayantara Gupta and Bing Zhang 



CO 



(J 
fat) 
M 
0) 

M 

0) 

M 

o 



io- 6 - 



io- 10 h 



10 



-14 _ 




z=0.1 
e e = 0.4,e B = 0.2 



i ■ i ■ i 



r = 800 
i.i.i 



10 4 10 8 10 12 
Observed Photon Energy (eV) 



10 



16 




-p 
Q 

'■p 

ft 

O 



3 



10 



10" 



-1 _ 




10" 

Observed Photon Energy 

Figure 5. An energetic short GRB. (a) The spectra of various components. Parameters: e e = 0.4, eg = 0.2, e p = 0.4, t v = 0.001s, V = 800, 
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the better targets for VERITAS and other TeV detectors. 
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